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A general approach allowing to construct the magneti- 
zation distributions of arbitrary topological charge in small 
exchange-dominated cylindrical ferromagnetic particles is pre- 
sented. The exchange energy functional is minimized by these 
distributions exactly. The magnetostatic energy is accounted 
partially, so that it facilitates a choice between the topologi- 
cally equivalent exchange-only solutions. The resulting mag- 
netization distributions can be easily generalized to a variety 
of non-circular cylindrical shapes by means of a conformal 
transformation. As an example a magnetic structures of a 
thin circular ferromagnetic cylinder both with centered and 
displaced magnetic vortex and of a finite Bloch domain wall 
in an elongated particle is given. 
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The magnetic structures of small ferromagnetic el- 
ements recently attained great experimental interest 
due to their possible applications in the integrated 
semiconductor-magnetic circuits and the magnetic ran- 
dom access memory (MRAM). From the point of view 
of these applications particularly the thin film elements 
(e.g. particles having cylindrical shape) are more inter- 
esting than others because they are more closely compat- 
ible with the current production technology. Currently, 
a lot of attention was devoted to modelling the mag- 
netic structure of these elements (see also references 
therein) with two main systematic approaches: the usage 
of guessed trial functions and the numerical solution of 
the Landau-Lifshitz equations. In this work an approx- 
imate (in accounting for magnetostatic interaction) but 
analytic approach for solving this task for a variety of 
element shapes is presented. 

Mathematically, the problem of finding the static 
distribution of the magnetization vector M(R) (R = 
{X, Y, Z} denotes the radius vector of a point in space) 
within an isotropic ferromagnetic particle in zero external 
magnetic field can be formulated as finding the minimum 
of the following energy functional of m{R) = M(R)/Ms- 




(1) 

with the constraint = 1 for R € V, where V 

denotes the volume of the particle, Le — y/C/Mg is 
the exchange length, C is the exchange constant of the 
material, M$ is the saturation magnetization of the ma- 
terial, V = {d/dX,d/dY,d/dZ} is the gradient opera- 
tor, hrj[m] = Hrj[m]/Ms is the demagnetizing field (a 



functional of m) created by the magnetization distribu- 
tion. The demagnetizing field can be expressed with the 
help of Maxwell equations through its scalar potential 
Ho = \U(K), which is, in turn, a solution of the Poisson 
equation V 2 t7 = V- M with the requirement (due to the 
finite size of the particle) that both \R\U and |.R| 2 |V£/| 
are finite as \R\ — ► oo, B. The main difficulty in treating 
this problem analytically arises from the fact that the 
dipolar interaction is non-local (as expressed by a func- 
tional dependence of ho on rh) , and the Euler-Lagrange 
condition for the extremum of the functional (|l|) yelds 
a system of non-linear three-dimensional integral partial 
differential equations. 

The simplified variant of the problem when the de- 
magnetizing field is assumed to be zero ho = (so 
that the distribution is defined by a system of partial 
differential equations, but still non-linear ones) and V 
is taken as infinite two-dimensional space with the con- 
straint m(i?) = mo at \R\ — > oo was solved by Belavin 
and Polyakov [0J|]. Let us briefly rederive this solution 
because the procedure will be relevant further in this pa- 
per. Noting that the constraint |m| = 1 leaves only two 
independent components of the magnetization vector, 
and parametrizing the magnetization vector by a com- 
plex function w(z,z) of a complex variable z = X + iY, 
~z = X — iY , i = V— T so that m x + im y — 2w/(l + uuw) 
and m z = (1 — ww)/(l + ww), the exchange energy func- 
tional can be written as 

2 f 4 / dw dw dw dw\ 2 
e = L E L z J v{1 + w _ )2 ^—— + ——j d R 2 , (2) 

where d/dz = (d/dX - td/dY)/2, d/&z = {d/dX + 
id/dY) jl and i?2 is a two-dimensional radius vector. The 
system of Euler equations for the extremum of this func- 
tional noting that the variables z and ~z are independent 
is 

d / dw\ 2w dw dw 
dz \dz J 1 + ww dz dz 

It is immediately clear that any analytic function w 
of the complex variable z satisfies this equation. Be- 
cause the Cauchy-Riemann conditions for the analyticity 
dw/dz = 0, if satisfied, turn both sides of (j^) to zero. 
The Belavin-Polyakov (BP) soliton corresponds to the 
representation of w(z) as a rational function of z. Zeroes 
of w(z) correspond to the centers of solitons (at these 
points m z — 1), at the infinity \z\ — > oo the BP soliton 
has m z = — 1. 
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Let us now turn to small particles in the form of cylin- 
ders, choosing the coordinate system in such a way that 
the axis Z is perpendicular to the faces of the cylinder. 
As long as the thickness of these cylinders is much smaller 
than the exchange length it is reasonable to assume that 
the magnetization distribution is uniform along the thick- 
ness (Z axis). This assumption immediately reduces the 
dimensionality of the problem and allows us to use the 
m(r) parametrization in terms of w(z,z). 

A special case of the problem ([l]) taking into account 
the demagetizing field was solved by Usov || for the par- 
ticle shaped as a circular cylinder and is also recently con- 
firmed experimentally pT0[ ] . This solution has the topol- 
ogy of a vortex, see Fig. . Its interesting property is that 
only a part of the magnetization distribution is a BP soli- 
ton (Usov calls this part the vortex core), the other part 
(represented in cylindrical coordinate system) as to,- = 1 
is not. The component of magnetization normal to the 
cylinder face m z changes from m z — 1 at the center of the 
magnetic vortex to m z — at the vortex core boundary, 
and continues to stay zero in the non-core region. 

The non-core part of the solution by Usov can be easily 
understood on the basis of the magnetic pole avoidance 
principle which is due to the fact that the energy 
of the demagnetizing field hu is always positive, and, 
since this field is created by magnetic poles (whose den- 
sity is V • M), the system tries to avoid having poles. 
At a boundary of a finite magnetic particle where \m\ 
abruptly changes from 1 to (no magnetic material) the 
density of magnetic poles is proportional to the compo- 
nent of the magnetization vector normal to the boundary. 
Thus, compared to the BP solution in the whole circu- 
lar cylinder the solution by Usov has less magnetic poles 
on cylinder faces and is more favorable from the point of 
view of magnetostatic energy. Unfortunately, the proce- 
dure employed by Usov can hardly be generalized to other 
particle shapes because it heavily relies on the cylindrical 
coordinate system. 

If we look back at equation (^|) we can easily suspect 
that it may have more solutions besides the BP solitons. 
Although it was shown that the BP solitons are the only 
solutions with the finite total exchange energy in the infi- 
nite space, this requirement becomes irrelevant as long as 
we consider finite particles. Another family of solutions 
of (0) was first found by Gross u% and is called merons: 



w S m(z,z) 



\f(z)\ < 1 



m/jmm \m\ > 1 



(4) 



they are expressed as wm(z, z) = f{z)/\J f(z)f(z), where 
f(z) is again an arbitrary analytic function. These so- 
lutions have to z = everywhere and do not produce 
magnetic poles on the cylinder faces. However, they di- 
verge at the points where f(z) — so that the exchange 
energy density is infinite at these points, therefore, the 
pure meron solutions can not be realized in the particles. 

We may note, however, that on the contour \f(z)\ = 1 
the meron continuously joins the BP soliton w(z) = f(z). 
Thus, the function 



is also a solution of (j3|) for any analytical function of com- 
plex variable f(z), is regular and continuous everywhere 
and has much less magnetic poles then the BP solitons. 
This solution for the magnetization distributions of small 
magnetic particles is the main result of this paper. By 
choosing the appropriate f(z) many realistic magnetiza- 
tion distributions can be obtained, some of the examples 
will be given later in the text. 

The solution (|^), generally, should have lower magne- 
tostatic energy then the corresponding BP soliton w(z) — 
f(z) because of less magnetic charges on the faces of the 
cylinder. However, such an advantage may turn out as 
unimportant for cylinders with a very small area of the 
region f(z) > 1. Unfortunately, the magnetostatic en- 
ergy is strongly shape dependent due to its long-range 
character and the estimates of this critical particle size 
need to be done separately in each case. Once the en- 
ergy of the solution (Q) is lower, the BP soliton is unstable 
with respect to the continuous transformation to (Q) by 
tilting spins in the region \f(z)\ > 1 towards the plane 
of the cylinder. Such a transformation does not have an 
associated energy barrier because wsm(z) and f(z) have 
the same topology. 

Another important fact, easily checked by direct sub- 
stitution, is that for any analytic function o(z), the func- 
tion w(o(z), o(z)) solves the equation (|J) if w(z, ~z) solves 
it. This property allows us to generalize the solution (^) 
using conformal transformations to particle shapes where 
the function f(z) is hard to choose on the basis of sym- 
metry argument alone. 

As an illustration consider a particle shaped as a circu- 
lar cylinder. From the symmetry it is natural to choose 
the BP solution in the form of a vortex centered in the 
cylinder (chosen as the origin of the coordinate system), 
which is f(z) = i{z/c)\ v \, where c is an arbitrary real 
number defining the scale and v is the integer topologi- 
cal charge of the vortex. In the case v = 1 this solution 
is equivalent to the one obtained and investigated both 
analytically and numerically by Usov and shown in 

Fig. . The definite value of c for a particular material and 
the geometry of the particle can be found by minimizing 
the total energy ([!]) with respect to c similarly as it is 
done in |Q. 

To displace the vortex the conformal transformation 
of the unit circle to the unit circle can be used, which is 
expressed by the rational function o(z) — (z— a)/(l — az), 
where a = X c + iY c is the new coordinate of the vortex 
center, \a\ < 1. The resulting configuration (|j) with 
f(z) — i(z — a)/(l — az)/c is shown in Fig. . It can be seen 
from the picture that this configuration corresponds well 
(at least qualitatively) to the application of an external 
field in the direction of Y axis. Also the case of \a\ — > 
1 corresponds to the case of the uniformly magnetized 
particle. The complete analysis of the energies (both 
Zeeman and the magnetostatic) and the stability of such 
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FIG. 1. The structure of the vortex with v = 1 and c = 1/2 
centered in the cylindrical particle. The outer solid circle is 
the contour of the particle, the inner circle is the boundary of 
the vortex core. 
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FIG. 2. The structure of the vortex with v = 1 and c = 1/2 
displaced by 1/2 with respect to the particle center. Note that 
the vortex core has the circular boundary but is otherwise 
asymmetrical. 



a displaced vortex in an external magnetic field will be 
published elsewhere fl3|| . 

Another example is a multi- vortex solution of an arbi- 
trary topological charge, which can be built by choosing 
f(z) in (^) as a product of powers (z — z.i) Ui , where Zi 
and Vi are the coordinate of the center and the topo- 
logical charge of z-th vortex, respectively. In particular, 
when there is an infinite number of Vi = 1 vortices sit- 
ting on a straignt line (say < z < 1) the resulting mag- 
netization distribution f(z) = zexp(/ log(z — x) dx) = 
(z — l)(z/(z — l)) 2 /e can be classified as a finite Bloch 
domain wall with ends at points and 1. The curves 
tangential to the vector rh are of elongated (but not ex- 
actly elliptical) shape. A similar distribution was indeed 
observed in elongated magnetic dots [[uj and interpreted 
as a two-domain configuration. 

To summarize, a general approach allowing to build 
the magnetization distributions with an arbitrary topo- 
logical charge in small exchange-dominated cylindrical 
magnetic particles of variety of shapes is presented. The 
produced magnetization distributions minimize the ex- 
change functional exactly and partially account for the 
magnetostatic interaction. There are two factors limiting 
the presented consideration to small particles only. The 
first is connected with the fact that the magnetization 
distribution is assumed to be uniform along the cylinder 
thickness. The second is that the exchange energy of the 
proposed solution within a circle of a given radius 
increases logarithmically with the radius. It means that 
the exchange energy diverges with increasing the particle 
size, thus other non-uniform configurations may become 
more energetically favorable in larger particles. 

The presented approach should coexist well with exist- 
ing numerical finite-element methods of micromagnetics 
as a way to specify the initial configuration with given 
topological properties, and to subject it to a final numer- 
ical refinement. 
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